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ABSTRACT 

We present sub-arcsecond imaging in [Ne V] of N49, the brightest optical SNR in the LMC. Between 
the "cool" optical and "hot" X-ray regimes, [Ne V] emission indicates intermediate temperatures for 
collisionally excited plasmas (2-6xl0 5 K), for which imaging has been extremely limited. We compare 
the flux in these images to the O VI measured spectroscopically by FUSE in individual apertures 
and find dereddened line ratios that are reasonably consistent with our predictions for intermediate 
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velocity shocks. The overall luminosity in [Ne V] for the entire remnant is 1.2 x 10' 
given the measured line ratios, implies an overall O VI luminosity of 1.5 x 10 38 erg s 
large radiative losses indicate that this material must have been shocked recently relative to the total 
lifetime of the remnant. We also explore the complex spatial structure. We find [Ne V] tracing the 
[O III] emission more closely than it does Ha, measure significant shifts (~0.1 pc) between the peaks 
of different emission lines, and find two orders of magnitude variations in the flux ratios for different 
filaments across the remnant. These properties as well as the general filamentary character of the 
optical emission suggest thermally unstable intermediate velocity shocks. 

Subject headings: supernova remnants: individual (N49) — ISM: structure — shock waves — ultraviolet: 
ISM 



1. INTRODUCTION 

The ultraviolet (UV) transitions of [Ne V] and O VI 
probe intermediate temperatures for collisionally ionized 
plasmas [(2-6) xlO 5 K] inaccessible in the optical or X- 
ray wavebands. This virtually unexplored temperature 
range has the most efficient radiative cooling and domi- 
nates the radiative losses in supernova remnants (SNRs). 
Of the potential diagnostics in this regime, the forbidden 
[Ne V] lines AA 3345.8, 3425.9A are the only bright emis- 
sion lines accessible from the ground. 

While the Far-Ultraviolet Spectroscopic Explorer 
(FUSE) has successfully provided emission line spec- 
tra from intermediate temperature SNR-shocks (Sankrit 
et al. 2001; Blair et al. 2002; Sankrit & Blair 2002; 
Sankrit et al. 2004; Korreck et al. 2004) it is not an imag- 
ing instrument. Study of the spatial distribution of this 
intermediate temperature material has been extremely 
limited. Low spatial resolution (~10') O VI images of the 
Cygnus Loop and the Vela SNR have been taken above 
the atmosphere with the High Resolution Emission Line 
Spectrometer (HIRES) (Rasmussen & Martin 1992) and 
SPEAR (Spectroscopy of Plasma Evolution from Astro- 
physical Radiation) (Nishikida et al. 2006) respectively. 
However, the only high spatial resolution (< 1") im- 
age to date of a SNR is the ground-based [Ne V] ob- 
servation of the Cygnus Loop (Szentgyorgyi et al. 2000). 
Ground based Fabry- Perot observations of 64 individual 
sites within the Cygnus Loop in [Ne V] that incorpo- 
rate spectral information from the tuneable filter were 
conducted by Sauvageot et al. (1999) but at much lower 
spatial resolution and hence lower sensitivity to sharp 
features. 

N49 is a good target for exploratory [Ne V] imagery 
for a number of reasons. It is the brightest optical SNR 
in the Large Magellanic Cloud (LMC), and hence has a 
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known distance, is likely to be detectable, and is very 
well-studied. Previous work has shown evidence of both 
radiative and hot X-ray emitting shocks. Shull et al. 
(1985) interpreted this in terms of the evolution of a 
SNR blast-wave into a complex circumstellar environ- 
ment sculpted by the progenitor star's radiation and 
winds propagating through an existing multi-phase in- 
terstellar medium. Given the presence of a much larger 
existing molecular cloud coincident with the south-east 
portion of N49 (Hughes et al. 1991), Vancura et al. 
(1992b) propose a simpler picture in which the temper- 
ature and velocity structure seen in N49 are the direct 
consequence of a Sedov-like blast-wave running into ex- 
isting large clouds. The filamentary structure of the op- 
tical image and the bands seen in the long-slit echelle 
data would both be the natural consequence of a "rip- 
pled sheet" geometry (Hester 1987) of alternating face-on 
and edge-on shocks as the blast wave envelopes the large 
clouds. Vancura et al. (1992b) derived a pre-shock den- 
sity of 20-940 cm -3 which they found to be inconsistent 
with Shull's model and state that their geometry better 
explains the structure of the echelle data and the fluxes 
seen. 

2. OBSERVATIONS, DATA REDUCTION AND ANALYSIS 

Observations of N49 were conducted on 2004 February 
16, at the Clay Telescope, one of the twin 6.5m Magel- 
lan telescopes at Las Campanas Observatory, Chile (-29° 
latitude) using the Raymond and Beverly Sackler Mag- 
ellan Instant Camera (MagIC). The CCD is a grade 
SITe 2048 x 2048 pixel array cooled to a detector temper- 
ature of -120° C. The 24/xm pixels each subtend 0.069" 
and the seeing varied between 0."6 and 0."8. At 50 kpc 
to the LMC each pixel corresponds to 0.0167 pc. We 
installed custom [Ne V], Ha, [O III], and associated con- 
tinuum filters manufactured by Omega Optical Inc. for 
our observations (see Table 1). The [Ne V] filter was nar- 



2 



Rakowski, Raymond, & Szentgyorgi 



TABLE 1 

Observations 



Filter 


Central A 


FWHM 


Exposure Time 


Airmasses 


Conversion 81 




(A) 


(A) 


00 




erg cm -2 counts -1 


[Ne V] 


3426 


15 


1200 


1.35, 1.38,1-42 


4.14E-15 


UV continuum 


3388 


35 


1200 


1.48,1.54,1.60 




Ha 


6563 


46 


100 


1.33,1.68,1.69,1.70 


6.18E-17 


red continuum 


7130 


58 


25 


1.34,1.72,1.72,1.72 




[O III] 


5007 


51 


50 


1.73,1.74,1.74, 


1.14E-16 








1.77,1.79,1.80,1.81 




green continuum 


5260 


59 


20 


1.82,1.84,1.86 





a Conversion from counts in median combined, continuum subtracted final images 



rower and the UV-continuum wider than those used at 
Whipple Observatory for the study of the Cygnus Loop 
(Szcntgyorgyi et al. 2000). 

SNR N49 was at an air mass less than 1.5 for the first 
half of the night before 4:00UT, during which time we 
concentrated on the [Ne V] observations. The Ha and 
[O III] images were taken primarily to facilitate compar- 
isons with [Ne V] . Table 1 lists the relevant details of the 
observations. Routine processing of the data was con- 
ducted in IRAF 1 making use of the MagIC package for 
the specifics of the CCDs. Flat-fielding was performed 
using twilight sky flats and cosmic-rays were eliminated 
by median-combining the separate frames. Given the 
extremely low surface brightness in [Ne V] and the UV 
continuum, high-frequency electrical noise in the CCD 
read-out which produced a "herringbone" pattern of 5 
to 10 counts could not be ignored. However, due to its 
slowly temporally varying periodicity it was possible to 
greatly reduce or even eliminate this signal relative to 
the diffuse target emission. 2 

Our data were calibrated with observations of the stan- 
dard star Hiltner 600 taken earlier in the night. The flux 
as a function of wavelength was taken from LCALIB in 
IRAF, in particular the flux of Hiltner 600 at 3426 A is 
4.35731 xl0~ 13 erg cm -2 s _1 A -1 at the top of the atmo- 
sphere. The extinction curves available in IRAF for the 
Cerro Tololo Inter-American Observatory (CTIO) were 
used, as they are a reasonable approximation to the ex- 
tinction at Las Campanas. The filter transmission curves 
had been measured at Omega Optical Inc., and these 
data are available at Las Campanas. 3 Over the [Ne V] 
and UV continuum bandpasses the extinction as a func- 
tion of zenith distance was well fitted by a third order 
polynomial, and the filter transmission by a Gaussian. 
For other filters simple linear or quadratic interpolation 
was preferred owing to features in the extinction curve 
or flat tops in the transmission as a function of wave- 
length. All other aspects of the throughput and detector 
efficiency were calibrated as a single conversion factor 
using Hiltner 600. The aggregate efficiency in [Ne V] 

1 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which arc operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation 

2 A scries of "PyRAF" scripts to identify and remove the pat- 
tern of electrical noise that combine the IDL and IRAF procedures 
written by S. Buries and M. Holman, is available on request to 
C. E. Rakowski. However owing to differences in IRAF and Python 
distributions and settings, "debugging" is likely to be required. 

3 The actual numbers were extracted from hard-copies of the 
transmission curves using a stand-alone version of "Dexter" from 
ADS. 



was sadly only 4.9%. We estimate that our photometric 
uncertainty is around 30% based on the mildly variable 
seeing and the small differences between the calibrated 
scaling factor for the continuum subtraction and the one 
which best subtracted the stars. Continuum subtracted 
images of N49 in Ha, [O III], and [Ne V] are displayed 
separately in Figure 1 and combined into a three-color 
image in Figure 2: Ha (red), [O III] (green), and [Ne V] 
{blue). 

The absolute coordinate system was found by compar- 
ison to the U.S. Naval Observatory (USNO) A2.0 cata- 
logue (Monet et al. 1998) using WCSTools (Mink 2002). 
At least 23 matched stars were identified for each filter 
with a mean radial offset of 0."367. However the error 
in the relative positions between the filters is a fraction 
of a pixel (0.2 - 0.75 pixel) since that is based on regis- 
tering the same stars in different frames using gcomap in 
IRAF. We note that the X-ray knot identified by Park 
et al. (2003) is ^5" south of the poorly subtracted star 
on the western boundary of Figure 2. Park et al. (2003) 
suggest that this knot is a fragment of ejecta from the 
SN explosion like those seen in the Vela SNR since it lies 
outside the apparent boundary of N49 and spectral fits 
indicate overabundances of O, Si and Fe relative to LMC 
values. There is nothing at the position of the knot in 
the [Ne V], Ha, or [O III] images that might be identified 
as cooler emission associated with the X-ray knot. 

3. DISCUSSION 
3.1. Comparison of [Ne V] and O VI 

The only information previously available on plasma 
at a few hundred thousand degrees in N49 comes from 
measurements of the O VI doublet with the Hopkins Ul- 
traviolet Telescope (HUT) (Vancura et al. 1992b) and 
FUSE (Blair et al. 2000; Sankrit et al. 2004) UV spec- 
trometers. These spectra provide important line profile 
information, but they are subject to uncertainties in red- 
dening and in the effects of resonant scattering by O VI 
ions, both within the emitting plasma (e.g. Long et al. 
1992; Raymond et al. 2003) and along the line of sight. 
There is also a question as to how much of the O VI 
emission arises from radiative shock waves in the 200 
km s _1 range, and how much comes from 300-400 km s _1 
non-radiative shocks that produce soft X-rays. Both this 
question and the uncertainty in the effects of radiative 
scattering within the emitting gas can be partly answered 
by examining the way the intensity is distributed within 
the apertures used for the O VI observation. 

The FUSE and HUT apertures are superposed on each 
of the emission line images of N49 in Figure 1. Table 2 
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Fig. 1.— Emission-line images of N49 Ho (left) [0 III] (middle) [Ne V] (right). HUT and FUSE apertures are overlaid. The HUT 
aperture is the large diagonal swath across the south-east while the 4"x 20" diagonal aperture is FUSE C05501. The 30" X 30" apertures 
are C05502, X00502, and C05503 from north to south. The images are linearly scaled over somewhat arbitrary ranges chosen to capture 
the full dynamic range of each image. The [Ne V] image was first binned by 4 pixels then smoothed with a 0.828" Gaussian kernel. Ho 
and [O III] were not binned or smoothed. 




Fig. 2. — Three-color image of N49 in Ho (red), [O III] (green), and [Ne V] (bJue). The [Ne V] image was first binned by 4 pixels then 
smoothed with a 1.38" Gaussian kernel. Ho and [O III] were not binned but were smoothed to 0.345". Regions from which the filament 
profiles were extracted are indicated. 
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TABLE 2 

Observed Intensities (10 -13 erg cm~ 2 s _1 ) 



Aperture 


[Nc V] 


Ha 


[O III] 


O VI a 


WFPC Ha a 


X00502 


2.19 


210. 


64.2 


9.9 


288. 


C05502 


1.69 


64.4 


28.0 


6.2 


99.0 


C05503 


1.39 


103. 


28.8 


5.1 


126. 


C05501 


0.34 


12.3 


6.60 


0.96 


20.0 


HUT 


1.30 


91.5 


33.6 


6.0 





a Sankrit et al. (2004) 



TABLE 3 

Dereddened Intensities (10 -13 erg cm~ 2 s _1 ) 



Aperture 


[Ne V] 


Ha 


[O III] 


O VI 


WFPC Ha 


X00502 


13.6 


549. 


221. 


2180. 


752. 


C05502 


10.5 


168. 


96.4 


1360. 


258. 


C05503 


8.66 


269. 


99.2 


1120. 


329. 


C05501 


2.13 


32.1 


22.7 


212. 


52.2 


HUT 


8.10 


239. 


115. 


1320. 





gives the measured intensities in the O VI doublet, along 
with the intensities in [Ne V] A3425, Ha, and [O III] 
A5007 measured by integrating over the apertures after 
masking out a few poorly subtracted stars. Table 2 also 
gives the Ha intensities reported by Sankrit et al. (2004) 
extracted from an archival HST image. 

In Table 3, we give the intensities corrected for a mean 
reddening E(B-V) = 0.37 as recommended by Vancura 
et al. (1992b). Vancura et al. (1992b) point out that the 
reddening is probably not uniform over N49, and that 
a scatter in the reddening value will introduce a correc- 
tion factor corresponding to a reduction in E(B-V). This 
is because some fraction of the area within the aperture 
sees significantly less reddening and contributes a dispro- 
portionate amount of flux. Sankrit et al. (2003) assume 
E(B-V) = 0.37 and that the scatter in E(B-V) is Gaus- 
sian, with a width of 0.14, which reduces the dereddening 
factor for O VI by a factor of 0.4. Sankrit et al. also es- 
timate that resonance line scattering within the remnant 
(Long et al. 1992; Raymond et al. 2003) attenuates the 
O VI flux by about a factor of 1.5, and that scattering in 
the halo of the LMC decreases the flux by another 20%. 
The resonance scattering partly compensates for the ef- 
fect of scatter in E(B-V), and the net result is to reduce 
the corrected flux of O VI by 28%. We include this factor 
in the fluxes given in Table 3. 

The first important comparison is the [Ne V]/0 VI 
ratio. We have computed a model with an updated ver- 
sion of the steady flow shock radiative code of Raymond 
(1979) and the LMC abundance set of Vancura et al. 
(1992a). The ratio of O VI to [Ne V] depends only weakly 
on shock speed in this range, and because the cooling is 
so heavily dominated by O and Ne at temperatures of 
3-6 xlO 5 K, only the relative abundances of O and Ne 
come into play. The Vancura et al. LMC abundance 
set gives O/Ne = 5.0. The major limitation of these 
models is their neglect of the thermal instabilities that 
affect shocks faster than about 150 km s _1 (e.g., Innes 
1992; Sutherland et al. 2003). However, the [Ne V] and 
O VI lines are formed at such similar temperatures that 
thermal instabilities should not change their ratio signif- 
icantly. The predicted intensity ratio of [Ne V] to O VI 



ranges from 0.0090 to 0.011 for shocks in the 200 to 300 
km s _1 range. The ratios derived from Table 3 range 
from 0.0062 to 0.010, in very good agreement with the 
predictions. This confirms that the reddening and reso- 
nance line scattering correction factors of Vancura et al. 
(1992b) are reliable. 

There is one major difference between the O VI and 
[Ne V] formation temperatures. The emissivity of O VI 
has a high temperature tail extending to about 2 x 10 6 
K, while [Ne V] does not. Hence enhanced O VI might 
indicate some emission from higher temperature shocks. 
In fact, the smallest [Ne V]/0 VI ratio is in FUSE aper- 
ture X00502 where the X-rays are particularly bright. 
However, the uncertainties in the measurements, the red- 
dening correction, the relative abundances of Ne and O, 
and the model calculations can easily account for the ob- 
served [Nc V]/0 VI ratios which lie about 20% below 
those expected from the radiative shock models. Thus 
the agreement of the observed flux ratio with the predic- 
tions indicates that the bulk of the O VI emission arises 
from 200-300 km s _1 radiative shocks as opposed to the 
non-radiative shocks that are responsible for producing 
the X-ray emission from N49. 

The second question is the degree of structure of the in- 
termediate temperature gas within the FUSE apertures. 
The [Ne V] image in Figure 1 shows that in each of the 
30" x 30" apertures most of the emission comes from 
filaments a few arcseconds across that cover perhaps 1 /3 
of the area within the aperture. For example, in the 
X00502 aperture the three brightest knots covering less 
than 2% of the aperture contribute 1/8 of the flux. This 
is what one would expect from the above result that the 
O VI arises from radiative shocks. When it is possible to 
obtain an actual O VI image, it should show less contrast 
than does the [Ne V] image, because scattering will se- 
lectively dim the bright, edge-on filaments in resonance 
lines (Cornett et al. 1992). 

The image of N49 gives a total [Ne V] luminosity of 
the remnant of 1.2 x 10 36 erg s _1 . From the [Ne V]/0 VI 
ratios, this implies L vi = 1-5 x 10 38 erg s _1 . This is 
8 times the X-ray luminosity measured by Long et al. 
(1981)(1.9 x 10 37 crg s" 1 from 0.15-4.5 keV). According 
to the models described above, the O VI luminosity is 
about 10% of the total energy dissipated by a shock in 
the 200-300 km s _1 range. Thus the luminosity of these 
intermediate velocity shocks is around 10 39 erg s _1 and 
in the 6,600 yr estimated lifetime of N49 (Park et al. 
2003) they would have radiated 2 x 10 50 erg, a signifi- 
cant fraction of the nominal energy of a supernova. Ac- 
cording to the shock speed distribution of Vancura et al. 
(1992a), slower shocks also produce a comparable lumi- 
nosity. However, if the SNR has recently reached the 
shell of a bubble produced by the progenitor (Shull et al. 
1985) or the dense cloud to the East and South (Hughes 
ct al. 1991) after passing rapidly through low density 
gas, the age of the shocks could be significantly less than 
6,600 yr. Indeed, Park et al. (2003) find a time-scale 
as short as 100 years based on the ionization state of 
the X-ray emitting gas in the SE. It is also worth noting 
that Lo vi is 2.5 times the infrared (IR) luminosity at- 
tributed to thermal emission from dust (Graham et al. 
1987). Williams et al. (2006) suggest that much of this 
IR emission comes from emission lines rather than dust, 
in which case it is produced by intermediate velocity or 
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slower shocks. 

For comparison, estimates of the O VI luminosity of 
the Cygnus Loop range from 1.2 x 10 36 erg s~ x based on 
a rocket observation (Rasmussen & Martin 1992) to 4.6 x 
10 36 erg s" 1 based on a Voyager Ultraviolet Spectrometer 
observation (Blair et al. 1991). These estimates would be 
reduced by a factor of 2 for the recent distance estimate 
of 540 pc (Blair et al. 2005). Thus the O VI luminosity 
of N49 is 50 to 100 times that of the Cygnus Loop. 

Rasmussen & Martin (1992) find that the blast wave 
can account for the O VI luminosity of the Cygnus Loop, 
while Vancura et al. (1993) find that the blast wave ac- 
counts for only 10 to 15% of the O VI photons. Ras- 
mussen & Martin (1992) point out the overall morpholog- 
ical resemblance between X-ray and O VI emission, while 
Danforth et al. (2001) discuss a region in the Cygnus 
Loop where intermediate velocity shocks clearly domi- 
nate. The [Ne V] data favor the interpretation of inter- 
mediate velocity shocks for N49. The blast wave temper- 
atures in N49 exceed 0.4 keV (Park et al. 2003), which 
is too high to produce O VI in the high temperature 
tail. There could in principle be substantial quantities 
of gas at 1 to 2 x 10 6 K whose X-ray emission would 
be overwhelmed by that of the hotter plasma. Indeed, 
the detection of [Fe XIV] (Dopita & Mathewson 1979) 
requires some gas at these temperatures. However, the 
[Ne V] to O VI ratio derived above shows that this hot 
gas does not dominate the O VI luminosity. 

3.2. Spatial distribution and flux ratios of [Ne V], 
[O III], and Ha 

Vancura et al. (1992b) noted the existence of a rim 
of [O III] extending outside of the Ha emission. They 
interpreted this as a photoionized precursor. Our images 
(as well as HST) show that the [O III] emission actually 
peaks outside of Ha. However, based on the fraction of 
[O III] emission in the narrow spike in the [O III] echelle 
profiles Vancura et al. (1992b) find that the precursor 
accounts for less than 20% of the [O III] emission. Thus 
the exterior [O III] filaments must be from the radiative 
shock itself. This suggests a physical separation between 
the [O III] and Ha in the post-shock gas. In Figure 2 we 
can identify many exterior and interior filaments where 
there appears to be a distinct separation between [O III] 
and Ha. 

Efficiently radiating intermediate velocity shocks are 
widely thought to be subject to thermal instabilities for 
v s > 150 km s _1 , but magnetic (Innes 1992) or cosmic- 
ray (Wagner et al. 2006) pressure can suppress or damp 
the instability. Large variations in the line fluxes and 
ratios and wide separations between the emitting re- 
gions of different ions may indicate thermally unstable 
shock conditions. In the non-steady-state radiative shock 
models of Innes (1992), that include the "cushioning" 
effect of magnetic pressure support on the collapsing 
thermally overstable filaments, a peak in [O III] emis- 
sion ahead of the Balmer line emission was indicative 
of the presence of secondary shocks. Under the par- 
ticular shock speed and ambient conditions they chose 
(Vs = 175 km s _1 , n = 1 cm~ 3 , _B = 3^iG) the sepa- 
rations between peaks are typically 4-7"x(D/l kpc) -1 
but can at times be as large as 20"x(D/l kpc) -1 . Innes 
(1992) also finds that the crossings of multiple secondary 
shocks result in relatively short-term time variations in 



both the line-intensities and line-ratios with line-ratio 
fluctuations of greater than a factor of 20. They further 
predict an [O III] photoionization precursor extending 
45"x(L>/l kpc) -1 with intensities 15-20% of the peak 
[O III] emission which can account for the echelle spec- 
tra. 

At the distance to the LMC the calculations of Innes 
(1992) translate into separations as large as 0.4". How- 
ever the densities required to explain the flux are at least 
a factor of 20 higher than modeled by Innes (1992) while 
the lower LMC abundances will increase the cooling time 
by a factor of 3. Thus we would expect the separations 
to be a factor of 6 smaller than the predictions of Innes 
(1992) unless there was significantly more pressure sup- 
port either from a higher magnetic field or cosmic-rays. 

To quantify the separations and flux variations, 18 re- 
gions were chosen along "straight" sections of filaments 
where relatively clean profiles could be extracted in each 
filter. These include the clearest examples of spatial sep- 
arations plus a few profiles to represent the more typi- 
cal cases where the peaks in all bands essentially coin- 
cide. The profiles, labeled "1" to "18" from east to west 
in Figure 2, were fit with Gaussians to determine the 
peak amplitude, position and full-width-half-maximum 
(FWHM) of the filament in each band above a flat back- 
ground level. Of these, 13 profiles were well-characterized 
by this simple model, while others with wider profiles or 
that simply reached a plateau were not. Selected profiles 
are shown in Figure 3 and plots of the relative distances 
between peaks, the ratios of the integrated fluxes, and 
the relative widths of the filaments are shown in Figure 
4. From this investigation we find the following: 

• Any given ratio of fluxes [Ne V]/ Ha, [Ne V]/ 
[O III], or [O III]/ Ha varies over 2 orders of magni- 
tude (as measured by the parameters of the Gaus- 
sian fits), with the highest [Ne V]/ Ha ratio occur- 
ing in region 18 at the far southwestern edge of the 
remnant. 

• In most of the regions the widths of the filament in 
all bands are similar, but in regions 1, 6, 7, 11, 12, 
and 18 the [Ne V] emission appears more extended. 
Region 11 (along the southeast bright portion) in 
particular shows a wider profile in [Ne V] extending 
ahead of the other two. 

• Measured separations between the peak flux in 
[O III] and Ha are as large as 0.15 pc (0.62") in the 
extreme cases of regions 3 and 17, and at least 11 
regions show separations larger than 0.05 pc (0.2"). 
In over half of the regions the peak in [Ne V] lies 
between the other two, closer to the peak in [O III]. 
The exceptions are regions 2 and 12 where [Ne V] 
and Ha essentially coincide, and 5, 6, 11 and 14 
where [Ne V] leads [O III] by as much as 0.04 pc 
(0.17"). 

The flux ratio variations are in line with the idea of 
different filaments being at different stages of thermal 
instability. In fact, we find that steady flow shock mod- 
els cannot explain [O III]/Ha>2/3 or [Ne V]/[0 III] > 
1/4 for LMC abundances, both of which are seen. We 
also expect that the [O III] and Ha emission from cooler, 
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Fig. 3. — Intensity profiles across selected filaments, Ha (diamonds), [O III] (pluses), and [Nc V] (squares). Six profiles 3,4,10,11,12, and 
18, were chosen to demonstrate the range of properties seen among the 18 extracted. The fainter [O III] and [Ne V] have been multiplied 
by the factors shown in each plot. All profiles were extracted from the original unsmoothed, unbinncd data and arc displayed in terms of 
the average counts pixel -1 versus the distance in pc to give a better sense of the statistics. Conversion factors from counts to flux and the 
exposure times are as given in Table 1. 
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compressed gas should be narrower than the [Ne V] fil- 
aments which trace the 2-6xl0 5 K shocked gas. We do 
see evidence for this in some places. The positions of 
peak brightness in [O III], Ha, and [Ne V] are often sig- 
nificantly offset from one another by distances of order 
0.1 pc similar to the predictions for a much lower den- 
sity, solar abundance plasma (Innes 1992). Thus if the 
flux ratios and separations in N49 are to be explained 
by thermal instability significant magnetic or cosmic-ray 
pressure support would be required to slow the progress 
of the instability. 

The filamentary nature itself is evocative of thermally 
unstable shocks such as the "ribbons of high density" in 
the simulations of Blondin et al. (1998). The structures 
in N49 appear similar in scale and curvature to the late 
times of their high ambient density simulation (Blondin 
et al. 1998, Figure 11). Quantitatively an actual fila- 
mentary structure of condensed ribbons is not required 
because the surface brightness contrast between the fila- 
ments and the diffuse emission (factors of 5-10 in [O III]) 
is still consistent with a rippled sheet geometry (Hester 
1987). However the similarities are highly suggestive. 

4. SUMMARY 

N49 exhibits an extremely bright complex structure in 
[Ne V] . The radiative losses alone imply a recent interac- 
tion with dense material relative to the age of the rem- 
nant. The ratios of [Ne V] to O VI emission, which both 
occur at 2-6xl0 5 K, are consistent with emission from 
200-300 km s _1 shocks, confirming the dereddening re- 
lation of Sankrit et al. (2003) and implying that most 
of the O VI emission arises from 200-300 km s _1 radia- 
tive shocks. For [Ne V], [O III], and Ha, which occur at 
decreasing temperatures post-shock, we find ~0.1 pc sep- 
arations between the peaks of their emission in individual 
filaments and a factor of 100 variation in their relative 
fluxes, including a number of places where the [O III] /Ha 
or [Ne V]/[0 III] flux ratios are higher than allowed by 
steady flow shock models at LMC abundances. These 
properties are consistent with the presence of thermal 
instabilities, but at the densities implied by the total 
flux, additional magnetic or cosmic-ray pressure would 
be required. 

Thermal instabilities have been almost universally 
predicted to occur in radiative shocks above v s > 
150 km s _1 , however observational evidence for them has 
been severely lacking. Density fluctuations in the pre- 
shock medium can mask or mimic the effects of thermal 
instability and it can be difficult to distinguish between 
different departures from steady-flow. For instance in 
the Cygnus Loop, high [O III] /Ha flux ratios are inter- 
preted as evidence for a truncated or incomplete recombi- 
nation zone behind a recent shock interaction (Raymond 
et al. 1980, 1988). In that particular region it seems un- 
likely that thermal instabilities have taken hold given the 
smooth variation of properties over such an extended thin 
structure (Raymond et al. 1988). In the slightly faster 
shocks of N49, on the other hand, we are interpreting 
the variability over the entire remnant as snapshots of 
different stages in thermally unstable shocks. In N49, 
the spatial separation and large variability are the keys 
to interpreting the non-steady-flow flux ratios. However, 
without additional magnetic or cosmic ray pressure these 
structures would have been unrcsolvable. 
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Fig. 4. — Gaussian line-profile parameters for the selected re- 
gions, (a) Flux ratios of [Ne V]/Ha vs. [O III]/Ha measured by 
the area under the Gaussian and dcrcddcncd. (b) Width in pc of 
[Ne V] vs. H«. (c) Separation in pc, between the peaks of emis- 
sion, [O III] to [Ne V] vs. [Ne V] to Ha. Overplottcd in (a) is 
the line [Ne V]/[0 III] = 0.1. Only regions for which a Gaussian 
line-profile captured the peak position and width were included in 
these plots. Errors shown are based on (a) 30% error in calibration 
and the square-root of the number of counts, (b) ±2 pixels, and 
(c) ±1 pixel. 
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